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Unstable particles in finite volume: 

The broken phase of the 4-d 0(4) non-linear cr-model 
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According to a proposal of Liischer it is possible to determine elastic scattering phases in infinite volume from 
the energy spectrum of two-particle states in a periodic box. We demonstrate the applicability of this method in 
the broken phase of the 4-dimensional 0(4) non-linear cr-model in a Monte-Carlo study on finite lattices. This 
non-perturbative approach also permits the study of unstable particles, the a-particle in our case. We observe 
the cr-resonance and extract its mass and its width. 
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1. INTRODUCTION 

Recently Liischer derived a relation between the 
energy spectrum of two-particle states in a fi- 
nite box with periodic boundary conditions and 
elastic scattering phase shifts defined in infinite 
volume 0. Since two-particle energy levels are 
calculable by Monte Carlo techniques, this rela- 
tion opens the possibility to extract phase shifts 
from numerical simulations on finite lattices. The 
scattering lengths have already been examined by 
means of the asymptotic large- volume behaviour 
of the two-particle energy levels ||]. Furthermore, 
Liischer's relation for two dimensions has recently 
been used to extract the scattering phase shifts of 
the 0(3) non-linear cr-model [|j and to study reso- 
nance scattering for two coupled Ising systems . 
As a next step towards the determination of reso- 
nance parameters in QCD we show the appli- 
cability of the relation in the 4-dimensional 0(4) 
non-linear cr-model. We investigate this model in 
its broken phase, where the cr-particle is unstable. 

In order to study two-particle scattering in fi- 
nite (spatial) volume Liischer starts with a non- 
relativistic model considering two identical mas- 
sive bosons in a periodic cubic box. The box size 
is assumed to be larger than the interaction range. 
He shows that the energy spectrum is determined 



by the scattering phases at the energy eigenval- 
ues. For states below the inelastic threshold this 
result is then carried over to massive relativistic 
quantum field theories by invoking an effective 
Schrodinger equation. 

Since this whole reasoning only works in the 
absence of massless (physical) particles, we have 
to equip the Goldstone bosons with a mass before 
we can apply Liischer's formulas. So we add an 
external source term to the lattice action of the 
4-dimensional 0(4) non- linear cr-model: 
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(1) 



The scalar field is represented as a four compo- 
nent vector of unit length: = 1. 

2. THEORETICAL BACKGROUND 

Consider a system of two Goldstone bosons ("pi- 
ous" ) of mass TO^ with zero total momentum in 
a cubic box of size . The two-particle states in 
the elastic region are characterised by centre-of- 
mass energies W or momenta k defined through 
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W = 2\Jml + k^, k = \k\, (2) 

with W and k in the ranges: 

2m^ < T4^ < 4TO7r -i^ 0<fc/m^<\/3. (3) 
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They are classified according to irreducible rep- 
resentations of the cubic group. Their discrete 
energy spectrum Wj, j = 0,1,2,..., is related 
to the scattering phase shifts di with angular mo- 
menta I which are allowed by the cubic symmetry 
of the states 

In the subspace of cubically invariant states 
only angular momenta / = 0, 4, 6, . . . contribute. 
Due to the short-range interaction it seems rea- 
sonable to neglect alH^^O. Then an energy value 
Wj belongs to the two-particle spectrum, if the 
corresponding momentum kj — \/ — 
is a solution of 

<5o(fc,)--0(^^) +J7r. (4) 
The continuous function (j){q) is given by 

tan (-0(g)) ^ \ , <^(0) = , (5) 

2oo(l,q ) 

with 2^00(179^) defined by analytic continuation 
of 

2oo(.s,g^) = -^ E {n'-qT'. (6) 
V 47r J^l, 

This result holds only for < g'^ < 9, which is 
fulfilled in our simulations, and provided finite 
volume polarisation effects and scaling violations 
are negligible |l| . 

In our Monte-Carlo investigation we calculate 
the momenta kj corresponding to the two-particle 
energy spectrum Wj for a given size L, and can 
then read off the scattering phase So{kj) from 
eq. (^. In order to scan the momentum depen- 
dence of the phase shift we have to vary the spa- 
tial extent L of the lattice. 

3. HOW TO EXTRACT THE TWO- 
PARTICLE SPECTRUM 

The two-particle states in our model are also clas- 
sified according to the remaining 0(3)-symmetry: 
They have "isospin" 0,1,2. Since we expect the 
CT-resonance to be in the isospin-0 channel, we re- 
strict ourselves to that case. The corresponding 
two-particle operators on a lattice of spatial ex- 



tent L are defined according to 
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a=l p 

where fi{p), i — 1,2,... is some basis of wave 
functions and <I>2 ^ is the spatial Fourier trans- 
form of the lattice field 

<l^^,,-L-3^$e,e2"-^^/^ (8) 

X 

We want to compare two different kinds of 
cubically invariant wave functions. First we take 
plane waves, 

fiip)=S^,p2 , i=l,2,... (9) 

Our second choice is motivated by the periodic 
singular solutions of the Helmholtz equation (see 

Mp)^ip^-qfr' , z-1,2,..., (10) 

with values qt — {L/2TT)y^{Wi/2)'^ — using as 
input the plane wave energy spectrum Wi of some 
test runs. These latter wave functions should 
have an enhanced projection especially at ener- 
gies close to the resonance mass. 

We also take into account the cr-field at zero 
momentum, since it has the same quantum num- 
bers as Oi{t) defined in eq. (|^): 

Oo{t) = ^^^^L-'Y.'^i,. (11) 

X 

The eigenvalues of the two-particle correlation 
function matrix 

Q,{t) = { {0,{t) -0,it + 1)) 0,(0) ) (12) 

with i,j = 0, 1, 2, . . . should decay proportional 
to exp{—Wit). Up to exponentially suppressed 
corrections this also holds if the set of states is 
truncated at some value i = r such that Wr is 
above the inelastic threshold. In this way we get 
a system of independent states inside the elastic 
region |^ . In our simulations this requirement is 
fulfilled for r < 6 depending on k and L. 
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Figure 1. Single particle energy values E{k) di- 
vided by the lattice and the continuum dispersion 
relation, respectively, as function of k"^ 

4. TECHNICAL DETAILS 

In order to generate our configurations we used 
the reflection cluster algorithm in its multi-cluster 
version generalised to the case of finite external 
source J. The Fourier transform is com- 
puted by a 3-dimensional FFT. In the case of the 
plane waves we also used the 2-cluster method 
PI t o check our results for the correlation matrix 
(O). The hopping parameter n — 0.315 and the 
external source J = 0.01 are tuned such that the 
cr-mass lies in the elastic region. 

The simulations are performed on cylindrical 
lattices ■ T. On the HLRZ-CRAY Y-MP8/832 
we can only run lattices of sizes 16^-32 and 20^-32. 
The larger lattices 24^ • 32, 32^ • 40 are simu- 
lated on the ''Landesvektorrechner" SNI S600/20 
in Aachen. In order to check the dependence on 
T we have performed runs on a lattice 243-60, but 
without getting significant changes of the results. 

Instead of calculating the eigenvalues of C(t) 
we follow refs. [^ji), choose some small to and di- 
agonalise C-^^^ito)CXt)C-^/^{to). After block- 
ing the data we determine the mean values and 
errors using the jackknife method. Finally the 



energy values are extracted by a 1-parameter fit 
to the mean eigenvalues according to the ansatz: 
exp{-{t-to)W,). 

Variations of the rank (r-|-l) of the correlation 
function matrix have only little influence on the 
results, while to must be chosen small in order 
to guarantee numerical stability of the eigenvalue 
calculation. 

Furthermore we have developed a method to 
determine the upper end t^ax of the fit range 
dynamically. The eigenvalues are sorted such that 
the corresponding eigenvectors at successive time 
slices are as parallel as possible. This should be 
equivalent to a sorting with respect to the size of 
the eigenvalues, but due to statistical fluctuations 
at some time slice t = tmax + 1 this fails, thus 
determining t^ax- 

5. NUMERICAL RESULTS 

In order to control scaling violations, we flrst ex- 
amine the single particle energies extracted from 
the exponential decay of the propagator 

( E E ^^,P^ ^-P ,* *p-.o ) - , (13) 

a=l p 

with k^ — {2'K/L)^i. Fig. |^ shows that the data 
are well approximated either by 

E{k) « ^Jml + P (l - 0.06(1) P ) (14) 

or by the lattice dispersion relation {ku = 2 sin-^ ) : 

E{k)^Eiat{k) = 2sisiiih{^\/ml + P) . (15) 

Under the assumption that polarisation effects 
are exponentially suppressed by a factor of 
cxp(— m^L) we get for the lowest state (fc = 0) 
the infinite volume mass = 0.229(1). 

Fig. H shows our preliminary results for the two- 
particle energy spectrum. The points connected 
by solid lines are the numerical results, bounded 
from above by the inelastic threshold Am-^. We 
get about the same numbers for both kinds of 
wave functions (^, (^0|). The lowest value repre- 
sents the bound ground state, which is not rele- 
vant for the discussion of scattering states. The 
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Figure 2. Two-particle energy spectrum 

dotted lines correspond to the free energy spec- 
trum: 



W 



= 2 




TieZ^, (16) 



with = 0, 1, 2, . . . We see the expected trend: 
the first and second level are coming close to each 
other forming the so called "avoided level cross- 
ing" , such that we get an impression where the 
resonance plateau (dashed line) might be. 

For each lattice extent L and energy level Wi 
in the elastic region the corresponding scatter- 
ing phase shift is calculated according to eq. (|^). 
Fig. 1^ summarises the results. The vertical dot- 
ted line represents the inelastic threshold (sec 
eq. (D). 

Finally we perform a fit to the data according 
to the Breit-Wigner resonance formula (solid line) 



tan I 5 — 



W — ma 

r./2 



(17) 



determining the mass and the width of the a- 
resonance: 



3.07(1) 



0.19(2) . (18) 
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Figure 3. Scattering phase shift So in the isospin 
channel 

6. DISCUSSION AND CONCLUSIONS 

The resonance mass m„ is about the same as the 
mass we have obtained from a propagator fit. Be- 
cause of triviality, perturbation theory should be 
applicable in this model, but the tree level ap- 
proximation gives about three times the numeri- 
cal value for Fo- in eq. (18). This deviation could 
at least partially be caused by lattice effects. For 
further details the reader is referred to [^. 
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